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Abstract. This paper examines some solutions for confluent and double-confluent 
Heun equations. In the first place, we review two Leaver's solutions in series of 
regular and irregular confluent hyper geometric functions for the confluent equation 
and introduce an additional expansion in series of irregular confluent hypergeometric 
functions. Then, we find the conditions under which one of these solutions can be 
written as a linear combination of the others. In the second place, by means of 
limiting procedures we generate solutions for the double-confluent equation as well 
as for special limits of both the confluent and double-confluent equations. Finally, 
we present problems which are ruled by each of these four equations and establish 
relations among Heun equations and quasi-exactly solvable problems. 
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^' 1. Introduction 

This paper is concerned with solutions of Heun equations and possible applica- 
tions of such solutions. We consider only the confluent and the double-confluent 
Heun equations (CHE and DCHE, respectively), and one limiting case of each of 

P^ \ these. The solutions for the CHE come directly from the differential equation, while 

the solutions for the other equations are obtained from the solutions of the CHE 
by limiting processes. Initially, we briefly discuss each of these equations and their 

£\j . connections; some more details are found in previous works [HI [151 [IBJ ES]. Then, 

we outline the main features of the solutions and the structure of the paper. 

C — ■ The CHE [TQl UJ [34] , also known as generalised spheroidal wave equation 

in the form used by Leaver [25] reads (u) ^ 0) 

d 2 U dU 

z(z -*o)^2+ ( B i + B ^^~ + i B 3 ~ 27 iu(z - zo) + lu 2 z(z - z )] U = 0, (1) 

where Bi, 7] and uj are constants and z = and z = zq are regular singular points 
with indicial exponents (0, 1 + B\/zq) and (0, 1 — B2 — B\/zq), respectively. At the 
irregular point z = 00 the behaviour of the solutions, obtained from the normal 
Thome solutions [25j [33] , is given by 

lim U(z) ~ e ±^ 2 W-(B 2 /2)_ ( 2 ) 

z— >oo 

The singularity parameter zq may take any value and, when zq = 0, the CHE 
gives the following DCHE with five parameters [25] 

z 2 ^ + (£?i + B 2 z) ^ + {Bs~ 2t]ujz + uo 2 z 2 ) U = 0, {B x ^ 0, u ^ 0) , (3) 

where now z = and z = 00 are both irregular singularities {B\ = and/or uj = 
are degenerate cases |15|). At z = 00 the behaviour is again given by Eq. (j2"j). while 
at z = the normal Thome solutions afford 

lim U(z) ~ 1, or lim U(z) ~ e Bl/z z 2 - B \ (4) 

z^0 z^0 



The CHE and the DCHE admit a limit which changes the nature of the irregular 
singularity at z = oo, keeping unaltered the other singular points. This limit is 
obtained by letting that [151 [16] 



uJ 



0, rj — > oo, such that 2-quj = —q, (Whittaker-Ince limit) 



(5) 



where q is a constant. It is called Whittaker-Ince limit because Whittaker and Ince 
have used a similar procedure to get the Mathieu equation (|12p from the Whittaker- 
Hill equation (JTTJ) [TO] [21] . The Whittaker-Ince limit of the CHE is 



z(z 



d 2 U dU 

^oj— + (Bi + B 2 z)— + [B 3 + q(z - z )j U = 0, (q ^ 0) 



dz 2 



Ji' 



(6) 



(if g = this equation can be transformed into a hyper geometric equation), while 
the Whittaker-Ince limit of the DCHE is 



dz 2 



(B 1 +B 2 z)^ + (B 3 + qz)U 



0, (q ± 0, £i ^ 0) 



(7) 



(if g = and/or Si = the equation degenerates into a confluent hypergeometric 
equation or simpler equations |15j). Eqs. ([6]) and ([7]) differ from the CHE and DCHE, 
respectively, by the behaviour of their solutions at the irregular point z = oo, which 
now is obtained from the subnormal Thome solutions [33] , namely, 



Im^^e 12 '^ 1 ^^/^ 



(8) 



in contrast with the behaviour of original equations (normal Thome solutions). Eq. 
([7j) also results when we take z$ = in Eq. ([6]). 

The preceding equations and their connections are summarised in the following 
diagram which is a modified version of a diagram given in [16]. The upper boxes 
display the CHE and the DCHE. The lower boxes show the Whittaker-Ince limits 
corresponding to the CHE and DCHE, respectively. 



z(z -20)—^ + (#1 + B 2 z)— + 
dz z dz 

[B 3 - 2tjuj(z - z ) + uj 2 z(z - z Q )] U = 



ZQ->0 



,d 2 U , n n ,dU 
■— T + (B 1 +B 2 z)—+ 
dz* dz 



(B 3 



2f]UJZ + lo 2 z 2 ) U 



(9) 
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and 7/ — > oo, such that 2-qu; 



$ 



^d 2 U rn „ ,dU 
z(z -z )tt + ( B i + B 2 z)—+ 
dz* dz 

[B 3 + q(z - zq)] U = 



z -+0 



,d 2 U 



dU 



z-— T + (B 1 + B 2 z)—+ 
dz* dz 



(B 3 + qz)U = 



(10) 



These connections among equations having different types of singularities become 
fully effective only when solutions of the CHE admit both the Leaver and the 
Whittaker-Ince limits. Counter-examples are provided by Hylleraas [20] and Jaffe's 
|22] solutions which admit none of these limits, as we can see by using Leaver's form 
for such solutions [25] . 

Now we introduce the Whittaker-Hill and the Mathieu equations which are par- 
ticular cases of both the CHE and the DCHE [TO]. A trigonometric (hyperbolic) 
form of the Whittaker-Hill equation (WHE) is [21] 



d 2 W 
du 2 



+ K* 



?? i 2 - (p + 1)£ cos(2ku) + -^ 2 cos{Anu) 



W = 0, (WHE). (11) 



If u is a real variable, this equation represents the usual WHE when k = 1 and the 
modified WHE when k = i. On the other hand, the Mathieu equation has the form 

m 

d id 

-7-j + a 2 [a — 2k 2 cos(2au)] w = 0, (Mathieu equation), (12) 

where a = 1 or a = i for the Mathieu or modified Mathieu equation, respectively. 
Some details about solutions for the WHE and Mathieu equation regarded as CHE 
or DCHE are given in Ref. [16]. The Mathieu equation is also a particular case of 
equation ([7]) as shown in section 3. Incidentally, the original Whittaker-Ince limit 
[T9| I2T] is obtained when £ — > 0, p — ► 00 so that p£ = 2fe 2 , k = a and 1? = a in the 
WHE. This gives the Mathieu equation. 

On the other hand, the solutions for the Heun equations, and in particular for 
CHE and DCHE, take one of the following forms 



00 00 N 

^Ja n f n (z) = ^2 a n fn(z), ^a n / n (z), ^a n /n(z), (13) 

n n=—oo n=0 n=0 

where the series coefficients a n satisfy three-term or higher order recurrence relations, 
fn(z) is a function of the independent variable z, N is a non-negative integer and the 
symbol ^ n holds for summation running from negative to positive infinity. These 
are called, respectively, two-sided infinite series, one-sided infinite series and finite 
series. The finite series are also known as quasi-polynomial solutions, quasi-algebraic 
solutions or Heun polynomials. 

Expansions in two-sided infinite series are necessary to assure the series conver- 
gence when there is no free constant in the Heun equations. Thus, all the parameters 
of the CHE and DCHE which rule the time-dependence of Klein-Gordon and Dirac 
test-fields in some Friedmannian spacetimes [6l [14] are determined from conditions 
imposed on the spatial part of the wave functions [351 ES]- Similarly, in the scat- 
tering problem of ions by a finite dipole [25] or by polarisable targets [15] all the 
parameters of the radial Schrodinger equation are known. 

When some parameters of the Heun equations assume special values, one-sided 
infinite series truncate on the right giving expansions in finite series. These Heun 
polynomials may be useful to get solutions for quasi-exactly solvable (QES) problems 
[371 EEl 1321 HD] • In effect, a problem is QES if admits solutions given by finite-series 
whose coefficients necessarily satisfy three-term or higher order recurrence relations 
[23| ; in contrast, a problem is exactly solvable if its solutions are given by (gen- 
eralised) hypergeometric functions. Actually, there are QES potentials for which 
the Schrodinger equation leads to the CHE and DCHE as well as to the general, 
biconfluent and triconfluent Heun equations (see Appendix A). 

Excepting possibly the Heun polynomials, in general the solutions for the Heun 
equations do not converge for the entire range of the independent variable. Then, it 
is necessary to consider two or more solutions converging over different domains and 
having the appropriate behaviours at the singular points. It is also necessary to take 
into account the transformation rules which generate new solutions from a known 
solution (these rules result from substitutions of variables which preserve the form 
of the Heun equations but modify their parameters). 

We will start with a set of three solutions for the CHE, represented by two-sided 
infinite series which have coefficients that satisfy three-term recurrence relations. 
These solutions admit both the Leaver and the Whittaker-Ince limits and so we can 
generate sets of solutions for all the equations discussed above. Solutions obtained 
from this set by means of transformation rules are given in sections 2.4 and 3.3. 



More precisely, in Sec. 2 we take two Leaver's solutions in series of regular 
and irregular confluent hypergeometric functions for CHE and introduce another 
expansion in series of irregular confluent hypergeometric functions - see the solutions 
given in Eqs. (J33a| and the recurrence relations (|33bp and (|33c[) . The expansion in 
series of regular functions converges for any z, whereas the two expansions in series 
of irregular functions converge for \z\ > \zq\. From the properties of the three-term 
recurrence relations and of the hypergeometric functions, we shall find conditions 
which permit to write one solution as a linear combination of the others in the 
region \z\ > \zq\. These conditions also assure that the series coefficients of the three 
solutions are proportional to each other. 

We advance that the convergence of the expansion in regular functions over the 
entire complex plane does not dispense with the expansions in series of irregular 
functions. In fact, if the expansion in regular functions is a linear combination of 
the others, its behaviour at z = oo must be given by a linear combination of the 
two behaviours given in equation ([2]), since each expansion in irregular functions 
corresponds to one of those behaviours. However, if for instance iujz is real, one of 
the exponentials exp (±iuz) goes to infinity. In this case, to get a solution finite 
when z — > oo we have to use one of the expansions in irregular functions. 

The solutions for the limits of the CHE, given in the above diagram, are obtained 
by the same procedure used in Ref. |16j . where a set having only two solutions 
in terms of one-sided series of confluent hypergeometric functions was considered. 
Thus, in Sec. 3 we find a set of three expansions in series of Bessel functions for 
Whittaker-Ince limit Q of the CHE, the three solutions possessing exactly the same 
series coefficients. In Sec. 4 we find that the solutions for the DCHE are again given 
by series of confluent hypergeometric, while the solutions for Whittaker-Ince limit of 
the DCHE are given by series of Bessel functions once more. 

In summary, in Sec. 2 we deal with the CHE, in Sec. 3 with the Whittaker-Ince 
limit © of the CHE, and in Sec. 4 with the DCHE and its limit ([7]). Sec. 5 presents 
some conclusions and points to problems which are ruled by the equations of the 
above diagram. Appendix A shows some relations among Heun equations and quasi- 
exactly solvable problems, while Appendix B gives some properties of the confluent 
hypergeometric functions. 

2. The Confluent Heun Equation (CHE) 

In this section we review the two Leaver solutions in series of regular and irregular 
confluent hypergeometric functions for the confluent Heun equation [25] and intro- 
duce the extra expansion in series of irregular confluent hypergeometric functions. 
The fact that the expansion in terms of regular functions converges for any z (the 
expansions in terms of irregular functions converge for \z\ > \zq\) distinguishes the 
present solutions from the Leaver expansions in series of Coulomb wave functions 
[25] since the latter converge only for \z\ > \zq\ (if the series is doubly infinite). 

Firstly, in Sec. 2.1, we recall some features of three-term recurrence relations for 
the series coefficients and supply properties of the confluent hypergeometric functions 
which enter the solutions for the CHE and DCHE. We also write down a Barber- 
Hasse solution in power series since this is important to obtain finite-series solutions 
as well as to cover the cases in which the expansions in series of regular confluent 
hypergeometric functions are not valid. 

In Sec. 2.2 we analyse the set constituted by the expansions in series of confluent 
hypergeometric functions. This is called the fundamental set of solutions because, 
by means of the transformation rules it originates other sets of solutions for the CHE 
and, by way of limiting processes, it affords sets of solutions for the other equations 



given in the schema of the first section. We find the conditions under which one 
solution of the fundamental set is given as a linear combination of the others. After 
that, we truncate the two-sided series from below in order to get one-sided series 
solutions as well. In Sec. 2.3, we discuss the convergence of new expansion in series 
of irregular confluent hyper geometric functions and, finally, in Sec. 2.4 we show how 
we can generate new sets of solutions by using the transformation rules for the CHE. 



2.1. General remarks and the Barber-Hasse expansions 

If the coefficients of one-sided series solutions are denoted by b n , then the three- 
term recurrence relations have the form 



ao^i + A)&o = 0, a n b n+ i + f3 n b n + <y n b 



'n-l 



(n > 1) 



(14) 



where a n , n and 7 n depend on the parameters of the differential equation. This 
system of homogeneous linear equations has nontrivial solutions for b n only if the 
determinant of the respective infinite tridiagonal matrix vanishes. This demands the 
presence of some arbitrary parameter in the differential equation. Equivalently, these 
recurrence relations imply a characteristic equation given by the infinite continued 
fraction 1 25 1 



A) 



«07l "172 "273 
01~ 02- 03- 



(15) 



which must be satisfied in order to assure the series convergence. 

If 7n = for some n = N + 1, where N is a non-negative integer, the one- 
sided series terminates at n = N and, consequently, gives a finite-series solution with 
< n < N |2l. Thus the recurrence relations can be written in the form 



/ 0o "o 

7i 01 «i 
72 02 



V o 



o \ 



a- 2 



In-i 0N-1 a N-l 
IN 0N J 



( b \ 
h 

b 2 



bN-i 
\b N ) 



(16) 



If the elements «j , 0i and 7, of the previous matrix are real and if 

aa i+1 >0, 0<i<N-l, 



(17) 



then all the N + 1 roots of its determinant are real and different [2] . This theorem 
is important to determine a part of the energy spectra in the case of quasi-exact 
potentials. On the other hand, if a n = for some n = N, the series begins at 
n = N + 1, but in this case one may set n = m + N + 1 and rename the series 
coefficients in order to obtain a series beginning at m = 0. 

Now, suppose that there is a second solution with coefficients c n satisfying 



«oci + A)Co = 0, a n c n+ i + n c n + 7„c n _i = (n > 1), 
where n is the same as in Eq. (|14p . Then, if 

a n %+i = O-nJn+i, 



;is) 



(19) 

it follows from Eq. (|15p that both solutions have the same characteristic equation if 
n takes the same values in both series: in these circumstances, b n and c n in general 



are proportional to each other. Such proportionality requires the same range for n 
in Eqs. (|14p and (|18p because there are cases in which the relation (|19p is formally 
satisfied, but one solution is given by a finite series while the other is given by an 
infinite series. In these cases one series breaks off on the right and the other on the 
left. 

The previous remarks can be extended to doubly infinite (or two-sided) series. 
These expansions present a parameter v which must be determined from a char- 
acteristic equation if there is no free parameter in the differential equation, or can 
be chosen at will if there is a free constant. The recurrence relations for the series 
coefficients b n now take the form 

a n b n+ ± + f3 n b n + 7„6„_i = 0, (-00 < n < 00) (20a) 

where a n , f3 n , 7 n and b n depend on the parameters of the differential equation as 
well as on v. These recurrence relations lead to the characteristic equation [25] 

a a-i 7o «-2 7-1 «-3 7-2 , «o 7i «i 72 a 2 73 , nrt ,s, 

ft = ^ITir- + FFF"' (20b) 

If there is a second doubly infinite series having the recurrence relation a n c n +i + 
flnCn + lnCn-i = such that the condition (fT9]) is fulfilled, then both solutions satisfy 
the same characteristic equation (I20b[) . The two series are really doubly infinite if 
neither the coefficients of 6 n +i and c n +i nor the coefficients of 6 n _i and c n _i vanish, 
since such conditions assure that the summation extends from negative to positive 
infinity in both solutions. This requires that 

a n , a n , 7 n and 7 n do not vanish for any n, (21) 

a requirement which imposes constraints on the parameters of differential equation 
and on the characteristic parameter v, in the case of two-sided series. These con- 
ditions will be useful for studying the sets of two-sided solutions of the CHE and 
DCHE. In addition, if we choose v such that a_i = a.~\ = 0, then the series are 
truncated on the left since the summation begins at n = 0. Thus, we have 

a_i = (5-i = =4> one-sided series with n > 0. (22) 

Other restrictions on v and on the parameters of the Heun equations come from 
the properties of the special functions used to construct the series solutions. Thus, 
let us consider expansions in series of regular and irregular confluent hypergeometric 
functions for the CHE, denoted by <f>(a,c;y) and ^(a,c;y) respectively. These are 
solutions of the confluent hypergeometric equation [12] 

dyi + { V) dy 



y^T2+( c -y)-n- a v = o, (23) 



where the parameters a and c will depend on summation index n, on the parameters 
of the Heun equations and also on the characteristic parameter v in the case of 
two-sided infinite series. In fact, the following four solutions for Eq. (|23p 



vl(y) = ®(a,c;y), ipl(y) = ^{a,c;y), 

(24) 
l pl(y) = e yy l - c <i>(l-a,2-c;-y), ^(y) = e^ y 1 '* *(1 - a, 2 - c; -y) 

are all of them defined and distinct only if c is not an integer [12J. Furthermore, if 

a, c and c — a are not integer , (25) 



then any two of the solutions (J24H form a fundamental system of solutions for con- 
fluent hypergeometric equation [12]. The formula 

¥(a, c; y) = -^—^— $(a, c; y) + ^r^ y X ~ c $(a - c + 1, 2 - c; y), (26) 
1 (a — c + 1) 1 (a) 

gives the analytic continuation of ^ in terms of <£. The expression of <£ in terms of 
^ is obtained from the previous one by using the relation T(z)T(l — z) = it / sv~n.(ir z) . 
One finds [12], 



e i7re(a-c)TV c \ 

$(a,c;y) = "—, ^/ fr(a,c;y) + e^ fr(c - a,c; -y), (e = ±1) (27) 

1 (c — a) 1 (a) 

where the plus or minus signs are to be taken throughout following the conventions 
er = 1 if (-l) = e in , e = -1 if (-1) = e~ i7T =► (-l) c = e inc , (-l) c = e~ i7TC 



(by setting c = 1/2 we see that e = 1 if we take \f— 1 = i, and e = — 1 if \f— T = — £). 
The relation (|26[) allows writing <I>(a,c;y) as a combination of a regular and an 
irregular confluent hypergeometric functions; analogously, Eq. (|27p gives \&(a,c;y) 
in terms of regular and irregular functions. 

The general form of the expansions in series of confluent hypergeometric functions 
for the CHE and DCHE is 

U(<p i n (y))=e- i « z J2<<Pn(y)> ( 28 ) 

n 

where 7r^ denotes the series coefficients. In the fundamental set of solutions (|33ap 
the argument of the hypergeometric functions ip l n is y = 2iuz, and the parameters 
are a = (B2/2) — irj and c = n + v + B 2 . We use only the functions ifn(y), ^niv) an d 
<Pn(y) because <£>^(y) would lead to a solution whose domain of convergence excludes 
all the singular points of the equation. As (fn(y) is discarded, the formula (p7|) is the 
only one necessary to connect the three hypergeometric functions. Consequently, the 
conditions (|25p are replaced by 

a, c and c — a are not zero or negative integers. (29) 



The present conditions combined with conditions (|2ip allow to use Eq. (|27p in order 
to write one solution of the CHE or DCHE as a linear combination of the others. 

Now we write down a Barber-Hasse solutions in series of (z — zq) for the CHE 
[3] \TU[ 125] . They converge for finite values of z. The first solution is 



[^arber^ = ju>» £ Q W ( z _ ^ Q z \ = finite ) ( 30a ) 

n=0 

where recurrence relations for the coefficients are given by {a_[ = 0) 



z (n + B 2 + |i j (n + 1) a^ + 



n (n + B2 — 1 + 2ia>2!o) 



+B 3 + iivzo [B 2 + ' 



aW + 2w rt + ii) + -^-l a^i = 0. (30b) 



i Jl) _|_ Oh,) ( n . -I- in -I - — 1 ^ «.W . 

This expansion will provide solutions to the CHE for the case in which the Leaver 
expansion in regular hypergeometric functions, given in the following, is not valid, and 
admits both the Leaver and the Whittaker-Ince limits. Furthermore, it gives finite- 
series solutions with < n < N when irj + B2/2 = —N. However, the Whittaker-Ince 



limit of this Barber-Hasse solution does not admit finite-series solutions since the 
coefficient of a n -\ is a constant: 2iuj (n + ir/ + B 2 /2 — l)—*q. 

We may generate a group containing 16 sets of solutions for the CHE by applying 
to the previous solution the transformation rules given in Sec. 2.4. These solutions 
for the CHE give solutions for the Whittaker-Hill equation (WHE). In effect, if U{z) 
symbolises the solutions for the CHE, the solutions W(u) for the WHE (jlip are 
obtained by writing |16] 



W{u) = U(z), 



z = cos (ra), (k = 1, i) 



(31a) 



where the parameters of U(z) are given in terms of the parameters of the WHE by 



zo 



1, B^ 



1 

~2 : 



Bo 



1, ^3 



1 



[(? + !)£-<?], 



iuj 



VX] 



p + 1 



. (31b) 



For the WHE, the Barber-Hasse solutions become even or odd solutions with respect 
to the change of the sign of u [16] . 



2.2. The fundamental set of solutions 

First we examine the two-sided infinite series solutions and then we obtain the 
one-sided series. It is worth advancing that the solutions in terms of regular confluent 
hypergeometric functions ^(a, c; y) are not valid when the first parameter a is zero or 
a negative integer. This statement is true only if the solutions are written so that the 
parameter a does not depend on the summation index n, since the Kummer relation 
()B.ip gives another representation in which the first parameter of <&(a, c; y) depends 
on n. 

Defining the function <&(a,b;y) by [25] 



£/ >, r ( c_ a ) a/ n r (c- a) 

<$>(a,c;y) = K — . <S>(a,c;y) - 



r(c) 

the fundamental set (1281) reads 



re 



a a(a + 1) 9 

l!c 2!c(c + 1) 



(32) 



r- ,.(:) = e -^ 2 ^(-l)"6« *(^-tJ7,n + i/ + B 2 ;2»a;«) , 



S2 

2 



E/°°(z) = e -k* ^(-l) n &«# f — - «y, n + 1/ + B 2 \ 2iujz\ 



(33a) 



J^ c { ^\-2wz) x - n - v - B - i v(l + ir)-^,2-n-v-B 2 ; -2iujz ) 

n. ^ 



2 ' 



where the recurrence relations for bh, and Cn are 



(n + i' + I ) ( n + 1/ + mj + -y ) b^f! + 



n + i / )(n + f + B 2 — 1 + 2iwzo) 
+B 3 + wzo (b 2 + —J b$ +2iuzo(n + v + B 2 + — -l\b < £_ 1 = Q. (33b) 



and 



(n + v + !)<£[ 1 + 



+B3 



(n + v)(n + v + B 2 - I + 2iujz ) + iuzo \B 2 -\ I 

V z o J 

4 1} + 2iujz (n + v + B 2 + — - 1 J ( n + v + ttj + y - 1 J cj^ = 0. 



(33c) 



Apart from a multiplicative constant, the coefficients bh and oh are connected by 

c h 1] =T(n + v + i V + f) b ( n ] , (34) 

provided that the argument of the gamma function is not zero or negative integer. 

The solutions U\(z) and U^°{z) have been taken from Eqs. (166) and (167) of 
Leaver's paper [25] with v replaced by i/-\-B2- On the other hand, by using a Kummer 
relation given in Eq. (jB.ip . £/f°(z) is rewritten as 



U 



co {z)=e wz Y J C i n ) ^(n + v + ir ] +^,n + u + B 2 - 1 -2iujz\ . (35) 



Then, it becomes clear that this solution can be obtained by substituting n + u for 
n in the one-sided solution given in Eq. (33a) of Ref. [16J and by allowing that 
the summation runs from negative to positive infinity. The solution U\ converges 
for any z [25J, while both C/f° and Uf° converge for \z\ > \zq\. From the fact that 
$(a, c; 0) = 1 and from Eq. (IB.2J) . it follows that 

lim U x {z) ~ 1, lim U?°(z) ~ e^V^, lim U^{z) ~ e ^ z -«?-^. (36) 

Thus, two different behaviours at z = 00 are included in the solutions belonging to 
the same set. Notice that the parameter v does not appear in Eq. (|2J) because this 
was obtained from solutions given by one-sided series (Thome solutions). 

Since we are dealing with three solutions for a second order linear differential 
equation, now we establish the conditions to get one of these as a linear combination 
of the others in a domain where the three solutions are valid. From the recurrence 
relations (|33bj) and (133c[) we find that the three series are really doubly infinite if 

v, v + B<i + -j- and v + ir\ + -^- are not integers, (37a) 

since under these conditions neither the coefficients of bn and Cn nor the coefficients 
of b n _ x and c„_ x vanish, that is to say, the series do not truncate on the left or on 
the right. The last condition also assures that the series coefficients are linked by 
Eq. (|34p and in turn this implies that there is a unique characteristic equation. If, 
in addition to conditions (|37a|) . the conditions Q29|) are also satisfied, that is, if 

n + v + B2, =!?■ — irj and n + v + irj + -£ are not zero or negative integers, 

then Eq. (|27p may be used to prove that any of the three solutions is a linear 
combination of the others in a region where the three solutions are valid (\z\ > \zq\). 
Since n runs from —00 to 00, the above conditions are equivalent to 

U + B2, v + irj-\ are not integers; irj is not zero or negative integer, (37b) 

which repeat one of the conditions (|37ap . 

In short, to express one solution in terms of the others, the three solutions must 
be given by two-sided series and the formula (|27p for analytic continuation of the 
hypergeometric functions must hold. These are the general conditions which may be 
applied to any set of solutions generated from the first set through the transformation 
rules of the CHE. In fact, they are equivalent to the conditions (|2ip and (|29p . 

If U\ is a superposition of U^° and £/f° in the common domain of convergence 
(| z J > |zo), then the behaviour of U\ when z — > 00 must be given by a combination of 
the behaviours of C/f° and t/f°. However, for certain problems, one of the expansions 
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in irregular functions may be inadequate when z — > oo and, consequently, U\ becomes 
inappropriate as well. 

Only the restriction on the values of (B2/2)—irj cannot be satisfied by a convenient 
choice of v. This restriction also arises if we consider the solution U\{z) by itself, 
disregarding its connection with the other solutions. In fact, if {B 2 /2) — ir/ = —m 
(m = 0, 1, 2, • • • ), the hyper geometric function <J>(a, c; y) which appears in U\ becomes 
a polynomial of degree m with respect to its argument [lj and, then, the summation 
from negative to positive infinity is meaningless. The solution 

00 



Uf(z) = e~ iuJZ Y, d^(z - z ) n , (\z\ = finite) (38a) 

n=0 

where recurrence relations for the coefficients are (d_{ = 0) 

z (n + B 2 + — J (n + 1)4+1 + n(n + B 2 -l- 2iujz ) 

+B 3 - iLoz ( B 2 + y- J dg> -2iu(n-ir)+y- 1 J 4-i = 0, (38b) 

takes the place of U% when {B 2 /2) — irj = —m. Notice that Uf was obtained from 
the Barber-Hasse expansion (|3Uap by substituting (— a>, — rj) for {00, rf). Furthermore, 
even if {B 2 /2) — irj = —m, by means of the transformation rules we can find two-sided 
series expansions in terms of regular confluent hyper geometric functions which hold 
for this case. 

Now we consider the one-sided series solutions. From the recurrence relations 
(|33bp and ()33c[) we see that for truncating the three solutions on the left at n = the 
only choice of v common to the three solutions is v = - see Eq. (|22p . We rewrite 
these one-sided series solutions as 

00 (—\\ n J^> ( n \ 

Ut(z) = z- lu)Z Y. T(n + B 2 ) ®{-J- i ^ n + B 2 ;2iuz\ , 

ra=0 

Uf°(z) = e~ iu)Z Y bg) * -^ - irj, n + B 2 - 2iujz , (39a) 

n=0 ^ 2 ' 

Uf°(z) = e luJZ ^ c$ * ( n + ir] + — , n + B 2 ; -2iuoz , 

n=0 V 2 ' 

where the recurrence relations are {b_[ = c_{ = 0J: 



(n + 1) (n + ir, + ^) b% + £«&« + 2iuz (n + B 2 + ^ - l) b^ = 0, 



(39b) 



and 



(n + ljcl 1 ' + /#><£> 



f-l)(» + * + f-^ 



+2io;zo ( n + £ 2 + — - 1 ) ( n + irj + -^ - 1 ) t£li = 0, (39c) 



where 



/#> = n( n + £ 2 - 1 + 2iu;z ) + B 3 + iuz (b 2 + f£ 



According to the previous subsection, if irj + (B 2 /2) = —/(/ = 0, 1, 2, •• • ), the series 
in U\ and £/f° break off on the right and these solutions reduce to Heun polynomials 
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(0 < n < I), while the solution U^° truncates on the left (n > I + 1). However, there 
is no need of considering these Heun polynomials since the Barber-Hasse solution 
(I30ap also suplies finite-series solutions with the same characteristic equation. 
For the solutions (J39aj the series are infinite, with < n < oo, if 



B 2 H — L and ir\ + -# are not zero or negative integers, 



-0 



(40a) 



as we see from the recurrence relations (|39bj) and ()39cp . The solutions can be con- 
nected by means of Eq. (p7|) if 



B 2 , -£■ — %r\ and irj + -£■ are not zero or negative integers, 



(40b) 



as we see from conditions (129p . Then, the relation ch = T[n -\- ir] + (B 2 /2)]bh is well 
defined and Heun polynomials are excluded from (J39aj. Under the above conditions, 
the coefficients of the Barber-Hasse solution (|30ap are proportional to c n , that is, 



,(D 



-,(D 



(z ) n T[n + B 2 + (B 1 /zo)Y 



(41) 



This gives another reason to consider the Barber-Hasse solutions in conjunction with 
the expansions in confluent hyper geometric functions. 

Finally, notice that there is another type of finite series which is not included in 
the Barber-Hasse expansions. It occurs when 



B 2 + ^ 



B, 



~<> 



-I, 



I — 1, 2, 0, • • • , 



(42) 



in which case the three series in (I39aj) terminate at n = I. These Heun polynomials 
are important because they are preserved by the Whittaker-Ince limit, in opposition 
to the Heun polynomials obtained from the Barber-Hasse expansions. 



2.3. Convergence of the third solution 

For one-sided infinite series the convergence of £7f° have already been established 
in Ref. [16]. Next we show that the two-sided infinite series converges in both 
directions, that is, when n — > 00 and when n — » — 00. First we write the solution as 



tjr\z) = e^ Y. 4V- B2 -"-^n(y), 



(43) 



where 



*„('/) - * ( l + in--^,2-n-v-B 2 ;y }. // 



-2iuz. 



To determine the convergence of the series, we have to find the ratios 
cgi *„+i(y) 



lim — ttt , lim 



& V *n-l(») 



C ( n ] ^n{y) 



,(1) 



(44) 



For this, in the first place we divide the recurrence relations ()33cp by rich and retain 
only the leading terms, that is, 



1 + 



v+1 



11 



"n+l 
J 1 ) 



+ 



n + 2v + B 2 - 1 + 2iujz + O 



n 



-\-2iujzq 



on / -1 

n + 2v + -B 2 + — + ir] - 2 + O I - 

2 z 



n 



c (1) 
c (1) 



0. 
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The minimal solutions for this equation are 



lim 



-n+1 



n— »+00 JJ-) 



lim 



„(l) 

"n— 1 



-2io;zn 



1 
n 



1 / Bi B 2 

1 + -M/-1 + — + -^ 

n \ zq 2 



B 2 + v-3 



n 



(45) 



On the other hand, from Eqs. (JB.4I) we get 



,62 \ 



- ( n + 1/ + irj + — j * n +i(y) + (n + z/ + B 2 -l- y)*„(y) + y^n-i(y) = 0. 



Dividing this by n^ n we obtain 



1 ( B 2 

H — f + in 4- — 
n V 2 



* 

* 7 



+ 



1 + 



u + B 2 -l-y 



n 



n *„ 



0. 



Then we can verify that 

lim — ^ — = 1 [ 1 + iri I , lim — ^ — 

n— »+oo W„ n \ 2 ) rc— >- ac w„ 



-(n + i/ + B 2 -l). (46) 



In fact there are other possibilities, but the preceding are the only ones compatible 
with Eqs. (B.5) and (B.6), respectively. Hence, since y = —2iujz, we find 



lim 



lim 

n— »— 00 



cj+x Vn+l(y) = £0 

c!2i y *«-i(y) 



i + i^ 2 + ^- 2 ) 

n V Zq J 



2 

1 + -. 
n 



(47) 



c { n ] Mv) 

Therefore, by the ratio test |41j the series in Eq. (|43p converges in the region \z\ > \zq. 



2.4. Other sets of solutions for the CHE 

In order to generate other solutions for the CHE we use the transformations rules 
resulting from substitutions of variables which leave the form of the CHE unaltered 
but change its parameters. Thus, if U(z) = U{B\, B 2 , B 3 ; zq,uj, 77; z) denotes one 
solution of the CHE in the Leaver form ([3]), we have the rules Ti, T 2 , T 3 and T4 
which operate as [16] 



T x U{z) = z 1+s i/*°£/(Ci, C 2 , C 3 ; zq,u, 77; z), 

T 2 U{z) = {z- zo) 1 - B2 - Bl / z °U(B u D 2 , D 3 ; z , u, 77; z), 

T 3 U(z) = U(B 1 ,B 2 ,B 3 ;z ,-u>, -77; z), 

TiU(z) = U(-B 1 - B 2 z ,B 2 ,B 3 + 2rju;zo; z ,-uj,t]; z - z) 



Zqt^O, 

\/z , 
\/z , 



(48a) 



where 
C x 
D 2 



-B 1 -2z , C 2 = 2 + B 2 + ^-, C 3 



Bo 



2B 1 



2gi 

Bi ( Bi 



B 3 +[l + 



B 2 + 



20 



±£±, D 3 = B 3 + ZL{^L + B 2 



1 



(48b) 



Applying these rules to the basic set (|33ap we may generate a group containing 16 
sets of solutions for the CHE. We will write only the subgroup obtained by using the 
rules T\ and T 2 in this order: 



Ti 



T-2 



Ti 



{U X ,U?,U?) A {U 2 ,U?,U?) A (U 3 ,U^,U^) A (U±,U?,Ur) (49) 
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where [Ui, [/f°, t/f°) denotes the first set of solutions (|33ap . 

In each of the following sets, the conditions for doubly infinite series are obtained 

(i) (i) 

by choosing v such that the coefficients of & n and c„_ 1 do not vanish in the recurrence 
relations. Besides this, one solution can be expressed as a linear combination of the 
others, if a, c and c — a are not negative integers, where a and c are respectively 
the first and the second parameters of the confluent hypergeometric functions: these 
conditions lead to some restrictions on the values of v as well as on the parameters 
of the CHE. Despite the notation, it is understood that the parameter v may be 
different in each set of solutions. 

The second set of solutions admits the Whittaker-Ince limit but does not admit 
the Leaver limit. It reads 

U 2 {z) = e ~ iuz z l+ ^y. 

„ \ z 2 zq J 

U$°(z) =e~ iujz z 1+ ^x 



Zq 2 Zq J 



U%°(z) =e iujz z 1+ *o x 



42) 



Efo\ I B-\ Bo „ 2B\ \ 

c£>V (n + v + 1 + in H + —,n + v + 2 + B 2 + ; -2iujz 

V ZQ 2 ZQ J 

(2) 

where the recurrence relations for c n are 

(n + v + l)c^ +l + (n + u)(n + u + l + 2iuz + B 2 + — - J 
+B 3 +(l + ^) (b 2 + *) + iuz (B 2 + *) 

\ Z J \ Zq J \ Zq J 

+2iujz In + v + B 2 + — - 1 J (n + u + in + — + ^ J c^_ x = 0. (50b) 

(2) 

The recurrence relations for b n are obtained from the previous ones by 

^ =T[n + u+l + in + ^ + ^-] b { n ] . (50c) 

The third set of solutions, which admits both limits, is 

U 3 (z) = e~ iujz (z - Zq) 1 ' 82 '^ z 1+ ^ x 

^(-l) n ^ 3 )$ (2-i V -^,n + u + 4-B 2 ; 2itvz) , 

U$°(z) = eT iwz {z - Zof~ B2 ~^ z 1+ ^ x 



B 2 A _ n . \ (51a) 



^(-l) n ^ 3) ^ ( 2 - in -,n + v + 4 - B 2 ; 2iioz\ 



I7f°(z) = e tuJZ (z - Zq) z « z + z o x 

B 2 



^c^ 3) ^ (n + u + 2 + in- — ^,ra + v + 4 - B 2 ; -2iu;z j 
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where the c^ 3) satisfy (CORRIGIR) 
(n + v + l)c% 
+ 



(n + v) (n + v + 3 - B 2 + 2iw,zo) + iwz ( 2 — J5 2 ) + B 3 + 2 - B 

+2iujz (n + u + l-B 2 -J (n + u + l + i7 1 --^-jc^l 1 =0 



r (3) 



(51b) 



(3) 

and the recurrence relations for 6„ follow from 

cl 3) =r[n + ^ + ir/ + 2-^] b ( n ] . (51c) 

The fourth set of solutions is 
U 4 (z) = e~ iuJZ (z - zq) 1 ' 3 ^^ x 

„ V z o l z o J 

x 

o X 

B x B 2 , , „ „ 2^1 „. \ (52a) 

UF(z) = e tuJZ {z - z o y~" 2 ~~ x 



U%>{z) = e~ iuJZ (z -zq) 1 Bz -o x 



„ \ Z 2 Z J 



a-Ba-U". 



V «4 4) * ( n + v + 1 + *»/ ---^,71 + ^ + 2-^2 -; -2iwz ) , 

„ V 2 2 z J 



where the recurrence relations for Cn are 



(4) 



(n + v + 1)^ + 



(n + v) [n + u + 1 — B 2 \- 2iujzo ) 

V zq J 






Bx f B x \ ( Bi 
+B 3 + — [B 2 + — -1)+ ioozo [2-B 2 

zo \ zo J V Zq 

+2iujzo ( n + v + 1 - B 2 - — ) ( n + v + irj - ^ - — ) cj 4 ^ = 0. (52b) 

V z J V 2 Zq J 

The b n and c n are connected by 

c { n ] =T[n + u + ir l +l-^-^}b { n ) . (52c) 

One-sided infinite series result when we take v = in the two-sided series solu- 
tions, since this restricts the summation to n > 0. 



3. Whittaker-Ince Limit of the Confluent Heun Equation 

In this section we show that for the Whittaker-Ince limit of the CHE, that is, for 
equation 

d 2 U dU 

z(z-zo)^ T + (B 1 + B 2 z)— + [B 3 + q(z-z )]U = 0, (q / 0), (53) 

the solutions of the CHE reduce to expansions in series of Bessel functions. The 
procedure is the same used in [16J for one-sided solutions, but we correct the follow- 
ing systematic error: the Bessel functions of the first kind J\ which appear in the 
solutions of Ref. [16] must be replaced by (— l) n J\. 
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In Sec. 3.1 we write the first set of solutions for Eq. (|53j) . We find that the 
three solutions have exactly the same series coefficients and no restriction must be 
imposed on the parameters of Eq. (153f) in order to write one solution as a linear 
combination of the others. In Sec. 3.2 we show how these solutions are obtained 
from the solutions of the CHE by using the Whittaker-Ince limit ([5]). In Sect. 3.3 
we use the transformations rules for Eq. (|53p to generate other sets of solutions. 
The limits of the Barber-Hasse solutions are given in [16J, but such expansions are 
unsuitable to get finite-series solutions. 

3.1. The first set of solutions 

The Whittaker-Ince limit of the fundamental set of solutions given in Eqs. (j33aH 
yields (i = 1, 2) 

U,{z) = ^(-l)"c«(^i)- (n+ ^ 2 - 1) J„ + , +B2 -i(2^), Vz, 

n (54a) 

U?(Z) = Y,^M\^T (n+U+B2 ~ l)H rX» + B^^), \Z\ > N, 
n 

where the limits of the recurrence relations (|33bp and (|33c[) are 

(n + v + l)c£k + [{n + u){n + v + B 2 - 1) + B 3 ] c« 

+qz [n + v + B 2 + (Bx/zq) - 1] c^ = 0. (54b) 

In these solutions J\(x) denotes Bessel functions of the first kind of order A, whereas 
H\ (x) and H^ (x) denote Hankel functions of first and second kind respectively. 
The solution U\{z) comes from the solution U\{z) of the CHE, while the solutions 
denoted by Uf (z) follow either from U 1 (z) or U\°° (z). This set admits the Leaver 
limit zq — ► 0. 

The Bessel functions which appear in the solutions are all independent since their 
Wronskians are 



W[J x (x),H^(x)] = 2i/(nx), W[j x (x),H®(x)] = -2i/(nx), 
W[H^\x),Hf\x)} =-U/(irx). 

Then, the relation 

J x ( X ) = ^[H^(x)+H^(x)], (55) 

can be used to write each solution as a linear combination of the others in a region 
where the three solutions are valid. 

On the other hand, for a fixed A the asymptotic behaviours of the Bessel functions 
as Ixl — ► oo are [13] 



J\( x ) ~ y ^ cos { x - \^ - 1 71 ") . I ar s A < ?r; 

h[ 1] (x) ~ \/2/(ttx) e ^'k^-i^, -7T < arg x < 2vr; ( 56 ) 

Hf\x) ~ \/2/(ttx) e -i{*-^-h), -2tt < arg x < vr. 



Thus, for the solutions C/{ we find 



lim U[ 1] {z) ~ e 2i ^ z*-%"%, lim u[ 2 \z) ~ e - 2i v^ z *-%~*. 



2-^00 2— >oo 
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The behaviour of U\ when z — » oo is a linear combination of these due to Eq. (|55p . 
Notice that, by setting 

w(u) = U(z), z = cos 2 (<7u), (<r = l,i), 

,2 (57) 

z = l, Si = -i, 5 2 = 1, 5 3 = V-f, ? = fc 2 > 

in Eq. (|53j) . we obtain the Mathieu equation (|12p . Then, from the solutions (|54ap 
and (|54bp we get the following even solutions for the Mathieu equation 



w\(u) = y (— 1)"<4^ [fccos(oix)] '' I/ J n+;/ (2fccos(a"u)), Vu, 



re 



«,«(«) = V(-l) n c«[fecos(^)] _n_ ^^(2A;cos(cJu)), |cos(au)| > 1, 



(58a) 



n 



where the coefficients c n satisfy 



n- 



+^+l)c« 1+ [(n + ^) 2 + (k 2 /2) - (a/4)] c« + fc 2 [n + u - (1/2)] c^ = 0. (58b) 



In this set of two-sided infinity series solutions, the first solution converges for any 
u, in contrast with the usual two-sided solutions for the Mathieu equation which 
converge, all of them, only for | cos(au)\ > 1 [T|l32]. 

One-sided infinite series are obtained by putting v = in the two-sided series 
solutions, what restricts the summation to n > 0. 

3.2. Derivation of the solutions 

To compute the Whittaker-Ince limits, first we rewrite the solutions of the CHE 
in a form convenient for using the formulas (|B.7p and (|B.8p . Thus, we rewrite the 
solutions (|33aP as (q = —2rjuj) 

Ui(z) = e-™^ r( { -V n f-. *(%-iTi,n + v + B 2 ^), (59a) 

^ T(n + v + B 2 ) \2 it]) 



qz 
2', — i 



(59b) 



fj?(z) = e lu>z J2 E n (-qz) l - n ~ u - B2 T Urj + n + u + ^jx 



* I 1 + it] - ^,2 - n - v - B 2 ; -^ ) , (59c) 

2 IT] 



where the new coefficients are defined by 



(-l) n 6n 1} „ . {iv) n+U C { n ] 



n — i L e 

n — r,ri ■ / 7-. in\i ' re 



r[l - n - v - irj - (B 2 /2)] ' T[iri + n + v + (B 2 /2)] ' 

Inserting these relations into Eqs. (I33bp and (|33cP we find 

(n + v + l)D n+1 +PJpD n 

+2icoz [n + v + irj + (B 2 /2) - l] [n + v + B 2 + (B 1 /z ) - l] D„_i = 0, 

and 

( j_ i u \ n + u + ir,+ (gg/g) ] ,«(!)„ 

(n + ^ + 1) : £„+i + /^ ; -E„ 

-2t?u;zo [n + v + B 2 + (5i/z ) - l] E„-i = 0, 
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where 

Pn ] = {n + v){n + v + B 2 -l + 2iu;z ) + B 3 + iujz (b 2 + ^\ . 
Thence, we find that the Whittaker-Ince limits of these relations are identical since 

limc^ = \\mD n = limE^, (u> — > 0, 77 — ► oo, 27/cj = —9). 

Denoting by c. n the above limits, we obtain the recurrence relations (|54b|) . 
Now, by letting (—if]) —* 00 and using Eq. (jB.7[) we find 

$ f -y -ir],n + v + B 2 ; — J ->T(n + v + B 2 )(-y/qz) " " 2 J n +u+B 2 -i( 2 V^)- 

Thus the limit of Eq. (|59aj) is the solution Ui(z) written in Eqs. (|54ap . On the other 
hand, to obtain the limit of the solution (|59bj) we define L\(z) by 



Lx(z) := T (1 - i77 - n - v - %) tt (^ - ir,, n + v + 5 2 ; g 

Then, when (—irj) — » 00 the relation (jB.8[) gives 

Li(z)->2(±iy/qz) " " 2 K n+u+ B 2 -i(±2i^/qz) = 

*(!-»—*> (VP) 1 " ,1 -^*ff£ >/+flB _ 1 (2VP) , 



Z7re 



(i) 
where in last step we have used the relations among the functions K\ and H\ given 

in Eq. (|B.10p . Inserting these into the limit of the solution (|59bj) and supressing 

multiplicative factors, we find the solutions U\ (z) given in (|54ap . The same solutions 

follow from the limit of the solution (I59ch . Notice that these are formal derivations 

which involve some tricks. However, we may check the resulting solutions by inserting 

them into Eq. (J53l) . 



3.3. Other sets of solutions for Whittaker-Ince Limit of the CHE 

If U(z) = U(Bi,B 2 ,B 3 ;zo,q;z) denotes one solution (or set of solutions) for 
Whittaker-Ince limit ([6]) of the CHE, then the rules 3F\, 2? 2 and 37% given by [15] 



3\U(z) = z x + B ^ z «U(C x ,C 2 , C 3 ; z , q; z), z ± 0, 

3 2 U(z) = (z-z ) 1 - B 3- B i/ z °U(B 1 , D 2 ,D 3 ;z ,q;z), zq^O, (60) 

3 r 3 U(z) = U(-B 1 - B 2 z , B 2 , B 3 - qz Q ; z , -q; z - z), \/z 

can generate a group having 8 solutions. The constants C\ and Di are defined in 
Eqs. (|48bl ). Next we use only 3\ and 3? 2 , following the sequence 

(u l7 U?) ♦ *♦ (y 2 , J7«) A (17a, U^) A (J7 4 , Uf) , (61) 

where (Ui,U[ ) denotes the set of solutions written in Eqs. (|54ap and ()54bp . The 
conditions for having doubly infinite series are obtained by choosing v such that the 
coefficients of c„ and (yl_i do not vanish in the recurrence relations, as in the case 



of the CHE. Besides this, according to Eq. (|55j) one solution can be expressed as a 
linear combination of the other solutions. We will write only the solutions in series 
of Bessel functions of the first kind, J\; the others are obtained by replacing J\ by 
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the first and second Hankel functions H^ (i = 1,2). The solutions for the Mathieu 



equation are derived by using Eq. (|57|) . 

The first solution of the second set, which does not admit the Leaver limit zq — ► 0, 

is (Vz) 



i+#i 



-n— V— 1— B2- 



u 2 (z) = z^Y.(- 1 ) nc n ) (Vqt) ^~ ' 



n+i/+l+B 2 



_2B X {2 y /qB), (62a) 



where the coefficients satisfy 



(n + v + l)^ + 



(n + v) ( n + z/ + B 2 + 1 + — - j + B3 



+ (l + ^-)(i? 2 + ^)] C ( 2 )+^o(n + ^ + i? 2 + ^-l)ci 2 2 1 = 0. (62b) 

For the Mathieu equation the corresponding solutions are even, as those of the first 
set, since 



W2(u) = cos (au) /(— l) n c„ [ A; cos (era)] 
with 



-n—v—1 



J n +u+i{2kcos(au)), Vu, 



(n + z, + l) C i 2 ] 1 + 



( n + z ,)( n + ^ + l) + _ + _^ 



.(2) 



+fc a [n + i/-(l/2)]c£> 1 = 0. 

In third set, which admits the Leaver limit zq — ► 0, the first solution is 



U$(z) = z z o (z - z ) z o x 



£(-i) n 4 3 ) (V5i) 



-n— 1/— 3+B2 



Jn+i/+3-B 2 (V^) , Vz 



(63a) 



associated with the recurrence relations 

(n + v + 1)0% + {(n + is)(n + v + 3-B 2 ) + B 3 + 2- B 2 ] c® 
+qz [n + u+l-B 2 - {B x /zq)] c^ = 0. 

For the Mathieu equation, this gives the odd solution: 



(63b) 



w 3 (u) = sin(2cra)£(-l) n c4 3) [A;cos(cra)] 

n 

where 



-n—v—2 



J n+u+2 (2kcos(cru)), \ra, 



k 2 a 



(n + u)(n + 1/ + 2) + — - - + 1 

2 4 



.(3) 



(n + ^+l)ci 3 j 1 + 
+fe 2 [n + v + (1/2)] cjfli = 0. 
The fourth set of solutions, which does not admit the Leaver limit, gives 



tt 1 \ 1 \i--B2- — 

V6\Z) = (z — zq) z o x 



D- 1 : 



M*) 



92) 



-n— v— I+B2" 



*o J 



n+v+l-B 2 



_r„_^l( 2 \/9 z ~)> Vz ' ( 64a ) 
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with the recurrence relations 



(n + u + l)c% + 

Bi ( B x 

+— [B 2 + — -l 



(n + v)[n + v + l- B 2 
(4) 



2Bi 

ZQ 



+ B 3 



<•;/ ! -! (/ : (1 ( a + v + 1 - B 2 - |^ J c^i = 0. (64b) 



For the Mathieu equation, we find odd solutions: 
W4(tt) = sin (cm) > ^(— l) n c4 ' [k cos(au)] 

n 
with the recurrence relations 

(n + i/ + 1)^! + 

+fc 2 [n + i/ + (l/2)]cj l 4 2 1 = 0. 



J n +i/+i(2£;cos(<7u)), Vu, 






,(4) 



4. The Double-Confluent Heun Equation (DCHE) 

In Sec. 4.1 we obtain sets of solutions for the double-confluent Heun equation Q 
by applying the Leaver limit (zq — > 0) to solutions of the CHE. The expansions in 
series of regular hypergeometric functions converge for any z, whereas the expansions 
in series of irregular functions converge for \z\ > 0. We also give the conditions to 
express one solution in terms of the others and obtain one-sided series solutions by 
truncating the two-sided series on the left. Finally, in Sec. 4.2 we write the solutions 
for Eq. flT]), that is, for the Whittaker-Ince limit of the DCHE. 

4.1. Solutions for the general case 



For zq — ► the solutions (|33aj) are not affected formally, but their recurrence 
relations change. We find 



E7i(s) =e-*"*J2(-l) n &W $ 



I 7 ( : | = e - iuJZ ^2(-l) n b£^ (^--i V ,n + v + B 2 ; 2iuzj , 



B2 

2 

2 



ir],n + v + i?2! 2iujz I 



(65a) 



C/f 5 ^) 



J^c^* (n + v + ii]+— -,n + v + B 2 ;-2iujzj 



Jl) 



The recurrence relations for the c« are 

(n + i/ + ljc^j + [(n + i/) (n + v + B 2 - 1) + iuB x + B 3 ] c^ 

+2iuB 1 [n + v + it) + (B 2 /2) - l] <^L X = °> 

wherefrom we get relations for fe„ by means of 

C { n ] =T(n + U + i7 ] +^)b i n ) . 



(65b) 



The behaviour given in Eq. (|36p for the CHE is also valid for the present solutions. 
These three solutions are given by doubly infinite series if 



v and v + irj + -£■ are not integers. 



(66a) 
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If, besides this, 

v + B 2 is not integer, and -£■ — ir\ is not zero or negative integer, (66b) 

then any of these solutions can be written as a linear combination of the others by 
using Eq. (|27p . Again, as in the CHE, only the restriction on (B 2 /2) — irj cannot be 
satisfied by an suitable choice of v. In addition, if (B 2 /2)—ir] = —m (m = 0, 1, 2, • • • ), 
the expansion U\ becomes meaningless. For this case, instead of U\, we can use 

Uf(z) = e-^J2 d n ) (4-) > (kl= finite) (67a) 

where the recurrence relations for the coefficients are (d_{ = 0) 

(n + l)4+i + [n (n + B 2 - 1) - iuB x + B 3 ] d$ 

-2itoB 1 [n-irj+ (B 2 /2) - l] d£l x = 0. (67b) 

This Uf, which reduces to a Heun polynomial when (B 2 /2) —ir] = —m, was obtained 
by letting z -> in Eq. ([38aj) . 

A second set follows from the solutions ()51a|) and the corresponding recurrence 
relations (I51bj) . Using the L 'Hospital rule we find that 



z 1+ ^{z - zo) 1 -* 2 -^ = z(z - zo) 1 ^ (l - ^) z » - t*- 3 *^'*, (z - 0) (68) 



Si , „ Si 

Z 

and, hence, we obtain the solutions 



f7 2 ( z ) = z 2-B 2e -ic 2 +^ Y^{-l) n b { ^ ( 2 - %n - — ,n + v + 4 - £ 2 ; 2iwz J , 



tf°°( z ) = z 2-B 2e -i^+^ yV-l) n i& 2) tf [2-ir]-—,n + v + 4-B 2 ; 2iujz) , 

n V 2 J (69a) 

C7 2 °°(z) = z 2 - B 2 e JLJZ+ ^x 



^c^ 2) ^ in + v + 2 + i7] -,n + u + 4- B 2 ;-2iujz) 



(2) 
where the coefficients c n obey 

(n + v + l)cjf +1 + [(n + 1/) (n + v + 3 - B 2 ) + B 3 + 2 - B 2 - iuBx)^ 
-liuBx [n + v + ir] + 1 - (£ 2 /2)] cjfli = 0. (69b) 

(2) 
The recurrence relations for the b n are derived from these via the relation 

cf) = r [n + v + ir] + 2 - (B 2 /2)] b®. 

Now the conditions to assure that the three solutions are given by doubly infinite 
series are 

v and v + it] — # cannot be integers, (70a) 

and if, besides this, 

v + B 2 is not integer, and irj + -^ is not zero or negative integer, (70b) 
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then Eq. (|27|) can be used to express one solution in terms of the others. If irj + 
(B2/2) = 2,3, • • • , instead of U 2 we can use the solution obtained from ([75]) by the 
substitution (77, <J)— ► (— 77, — u>). 

One-sided series solutions are obtained by setting v = in the two-sided solutions, 
in the same way as in Sect. 2.2. We write only the first set and the limit of the 
Barber-Hasse solution (|30ap . 

00 (—l) n r^ I B \ 

U 1 {z) = e~^ z Y, v( + B) ®[^--^n + B 2 -2iuz\, 

Uf°(z) = e~^ z ]T(-l) n 6« tf (^ - irj,n + B 2 ; 2iuz\ , (71a) 

n=0 ^ ' 

Uf°(z) =e !U2 ^cW$ (n + ir} + -£,n + B 2 )-2iuz) , 

n=0 ^ 2 / 

where the recurrence relations for 6„ ' and c„ are (&_{ = c_j = 0) 

(n + 1) [n + iri + (B 2 /2)] b% + [n (n + B 2 - 1) + iw-Bi + B 3 ] b$ 

+2iuB 1 6^! = 0, (71b) 

(n + ljcj^i + [n(n + B 2 -l)+ iuB x + .83] c^ 

+2iooB 1 [n + Z77 + (S 2 /2) - l] c^ = 0. (71c) 



The limit of the Barber-Hasse solution (|3Ua|) is 

U^ he *(z) = e i »*J2 c £ ) 14-) ■ (72) 

n=0 ^ *' 

Notice that the solution U^°(z) does not admit finite-series solutions, while the other 
solutions do. However, for finite-series it is sufficient to consider only Barber-Hasse 
solution. On the other hand, the condition in order that the three expansions in 
series of confluent hypergeometric functions are given by infinite series (n > 0) is 

irj + -^- is not zero or negative integer, (73) 

since in this case the series do not truncate on the left or on right. In addition, if 

B 2 and -^ — irj are not zero or negative integer, (74) 



then Eq. (|27l) can be used to connect the three one-sided infinite-series solutions. 
The Barber-Hasse solution corresponding to the second set is |16] 



00 / v n 

[/b-ber (z) = e iu z+{ B l/z ) z 2-B 2 £ c (2) I _ _f_ \ ^ (75) 

n=0 ^ X ' 

where the recurrence relations for c« are obtained by putting v = in Eq. (|69bp . 
The corresponding solutions in series of hypergeometric functions can be obtained 
by taking v = in the solutions (|69ap . 
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4.2. Solutions for the Whittaker-Ince limit of the DCHE 

For the Whittaker-Ince limit (|7|) of DCHE, namely, 

z2 l^ + (jBl + B2Z) 1J + {Bs + qz ^ U = °' (« / 0, ^ / 0), 

there is no finite-series solutions. For this reason, the equation is irrelevant for quasi- 
exactly solvable problems. However, its solutions may be important for studying 
the scattering of ions which induce dipole and quadrupole moments in polarisable 
targets [15]. Solutions for the above equation can be found by applying either the 
Whittaker-Ince limit to the solutions given in Sec. 4.1 or the Leaver limit to the 
solutions (|54a|) and (|63a|) . We use the latter approach because it is the easiest. 

In effect the Leaver limit zq — ► does not modify the solutions (|54ap , which again 
read 



n 

vi%) = ^- l ) n ^^T {n+v+B2 ~ 1)H( ^ + B,-A^)i \z\ > o, 



(76a) 



but changes their recurrence relations to 

+ v + 1)4^ + [(n + v){n + v + B 2 - 1) + B 3 ] c« + qBy^ = 0. (76b) 



m 



These three solutions are really doubly infinite if v is not integer. Furthermore, each 
solution can be expressed as a linear combination of the others by means of the 
relation (l55j) . 

On the other side, the only detail to get the Leaver limit of the solutions given 
in Eqs. (I63ap is the use the limit written in Eq. (I68p . The results are 

U 2 {z) = A 2 ^ £(-ir C ( 2 ) (^qz-y n - U - 3+B2 J n+U+ ^ B2 {2^), Vz; 

(77a) 
U?\z) = e^z^ E(- 1 ) n ^ ) (v / ^)" n_ ^ 3+B2 ^t + 3- B2 (2v^), \z\ > 

n 

where 
n + v + l)c% + [{n + v)(n + v + 3-B 2 ) + B 3 + 2- B 2 ] c^ - qB^^ = 0. (77b) 
One-sided series are obtained by taking v = once more. 

5. Conclusions 

In this section we collect some results we have found, mention other solutions 
which deserve further analysis, and indicate the possibility of applying the solutions 
of the four equations given in the first section to physical problems. 

We have started with a set of solutions for the confluent Heun equation (CHE) 
given by three series of confluent hyper geometric functions and, by means of Leaver 
and the Whittaker-Ince limits, we have derived sets of solutions for all the equations 
included in the diagram described in the first section. The Barber-Hasse expansions 
in power series also admit the both limits and can be used to provide solutions for 
the cases in which the solutions in hypergeometric functions are not valid, as well as 
to get finite-series solutions for the confluent and double-confluent Heun equations. 

In the fundamental set of two-sided solutions (Ui, C/J° ', Ui°), the first and the sec- 
ond solutions are, respectively, Leaver's expansions in series of regular and irregular 
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confluent hypergeometric functions for the CHE [25], while £7f° is the two-sided ver- 
sion of the expansion in series of irregular functions given in Ref. [16] . We have seen 
that, although the solution U\ is not valid if there is a certain constraint between two 
parameters of the equation, we can use the transformation rules to find a two-sided 
solution valid for that case. We have also established the conditions to express each 
of the three solutions as a linear combination of the others. 

The fact that the two-sided expansion U\ in series of regular functions converges 
for any z, and in particular in a neighbourhood of z = 0, distinguishes the present 
solutions from the Leaver expansions in series of Coulomb wave functions [25] , since 
the latter converge only for \z\ > \zo\. On the other hand, we have seen that 

lim Uf°(z) ~ e~ iulz z ir >-^, lim U?°{z) ~ e^-^J-"-^ 

where v is the characteristic parameter of the two-sided series. The behaviour of U\ 
when z — > oo is given by a combination of these. This explains why U\ is sometimes 
inappropriate when the variable z tends to infinity. In this event we need to use only 
one of the expansions in series of irregular functions. 

All the solutions for the double-confluent Heun equation (DCHE), obtained from 
solutions of the CHE through the Leaver limit {z$ — > 0), are also given by expansions 
in series of confluent hypergeometric functions. Then, the analysis used for the 
solutions of the CHE has been promptly adapted to the solutions of the DCHE, 
including the conditions which allow to write one solution of a fixed set as a linear 
combination of the others. 

We have also found one-sided infinite-series solutions for the CHE and DCHE, 
by truncating the two-side infinite series on the left. The one-sided solutions for the 
CHE can afford two different types of finite-series solutions. One of these may be 
ignored because they are also suplied by the Barber-Hasse expansions in power series. 
The other type of finite series is important because it is the only one preserved by 
the Whittaker-Ince limit of the CHE. 

The Whittakker-Ince limit ([5]) of the CHE and DCHE have generated expansions 
in series of Bessel functions of the first kind and in series of the first and second 
Hankel functions. In this case, each solution belonging to fixed set can be written as 
a linear combination of the others without the need of restrictions on the parameters 
of the differential equations. The Whittakker-Ince limit of the DCHE does not admit 
finite-series solutions. 

In order to have a full account of the solutions, we have considered the transfor- 
mations rules which result from substitutions of variables which transform a Heun 
equation into another form of itself (up to a change of parameters). For instance, 
only if we apply the transformation rules to solutions of the CHE, can we derive solu- 
tions for the Mathieu and Whittaker-Hill equations having all the desired properties 
of parity and periodicity. Similarly, if an expansion is not valid for certain values of 
the parameters, the transformation rules allow to find another solution valid for that 
case. 

Solutions in series of Coulomb wave functions |25| and in series of Gauss hyper- 
geometric functions [14} [28l [29] for the confluent Heun equation also admit both the 
Leaver and the Whittaker-Ince limits [144115] . However, there are several unexplored 
aspects concerning such expansions. Thus, it is necessary to consider other solutions 
of this type in order to retrieve the usual solutions of the ordinary spheroidal equa- 
tion [H [21 [32], which is a particular case of Eq. (pQ) with rj = [25] . This would 
amount new solutions for the equations of the diagram written in the first section. 
On the other side, by means of the Whittaker-Ince limit these solutions conduct to 
known solutions for the Mathieu equations [15]. Nevertheless, the possibility of find- 
ing solutions for the CHE which reproduce all the known solutions of the Mathieu 
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equations is an open issue. 

Now we refer to some problems which obey the CHE and its limiting cases. In 
the first place, the Schrodinger equation with inverse fourth and sixth-power poten- 
tials leads, respectively, to the DCHE ([3]) and its Whittaker-Ince limit ([7]). These 
potentials have appeared, for example, in the description of intermolecular forces [17] 
and in the scattering of ions by polarisable atoms [23] . 

In fact, the radial equation for the scattering of ions by neutral targets with 
induced dipole and quadrupole moments can be put in the form of the Whittaker- 
Ince limit of the DCHE [15]; for neutral or charged targets presenting only induced 
dipole moment, the equation may be reduced to the DCHE [3 [15]. In these problems, 
there is no arbitrary constant in the differential equations and so, in order to assure 
the series convergence, it is necessary to use two-sided infinite series solutions with 
a characteristic parameter v. 

In the second place, as already mentioned, the equations which rule the time- 
dependence of Klein-Gordon and Dirac test-fields in some nonflat Friedmannian 
spacetimes can be written as CHEs and DCHEs [U [14]. These equations also re- 
quire two-sided infinite series solutions. In addition, there are the CHEs for angular 
and radial Teukolsky equations in Kerr spacetimes [U [25[ [28l I29| . 

Finally, in Appendix A we have found that normal forms of the Heun equations 
are helpful to know if a given one-dimensional Schrodinger equation reduces to Heun 
equations. In this connection, we have seen that there are quasi-exactly solvable 
(QES) potentials representative of each of the five Heun equations, a fact which may 
become important if we have sufficient knowledge about the solutions of the equation. 

As noticed at the beginning, QES problems demand finite-series solutions which 
give only a part of the energy spectrum, but in the case of Heun equations there 
are also infinite-series solutions. Whether these afford the remaining part of the 
spectrum, this is a theme to be analysed. At last, we observe that the inverted QES 
potential (|A.7|) leads to the Whittaker-Ince limit of the confluent Heun equation. 
As far as we are aware, the only other example known for this equation (excepting 
Mathieu equations) comes from the separation of variables of the Laplace-Beltrami 
operator for scalar fields in an Eguchi-Hanson space [5[ I271I3T] . 

Appendix A. Heun Equations and Quasiexact Solvability 

We consider the Heun equations in their normal or Schrodinger form, 

" d 2 



dz 2+Q ^ 



y(z) = 0, (A.l) 



that is, in the form where there is no first-order derivative term |10| . The function 
Q(z) for the general Heun equation and its confluent cases is given in the following 
|34j . where it is understood that, in general, there exists some constraint among the 
parametres A, B, and so on. In each case we indicate the singular points of the 
equation. 

General Heun equation. Four regular singular points at z = 0, l,a, oo: 

«/ x A B CD E F , , . ,. . 

^ ) = 7 + — + 7^ + ^ + (73 T )I+(7^ F (^Oorl). (A.2) 

Confluent Heun equation (or generalised spheroidal wave equation). Two regular 
points at z = 0, 1 and one irregular at z = oo: 



B CD E 

z z — 1 z 2 (z — 1) 



Q(z)=A+- + - — - + — + - — t^. (A.3) 
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Double- confluent Heun equation. Two irregular points at z = 0, oo: 

n, S , B C D E /A ^ 

Q( z ) = A+- + - f + -3+- I . (A.4) 

Biconfluent Heun equation. One regular point at z = and one irregular point 

at z = oo: 

Q(z)=Az 2 + J Bz + C + - + 4- (A.5) 

z z z 

Triconfluent Heun equation. One irregular point at z = oo: 

Q(z) = Az A + 5z 3 + Cz 2 + Dz + £. (A.6) 

Other normal forms are given in the tables constructed by Lemieux and Bose 
[26] who, however, have not considered the triconfluent equation. These tables are 
helpful to recognise whether a given equation is of the Heun type. In particular, by 
using the Lemieux-Bose tables in addition to the normal forms written above, it is 
straightforward to establish relations among the Heun and the Schrodinger equations 
for some quasi-exactly solvable potentials. 

We find: (i) a triconfluent Heun equation for the quartic potential given in Eq. 
(5.34) of Gonzalez-Lopez, Kamran and Olver |18j : (ii) biconfluent Heun equations 
for the sextic potential V\z 6 + V22 4 + V32: 2 + V4 + V5/2; 2 given by Turbiner [38] 
and Ushveridze [40j . and for the potentials II, III and VIII given in Turbiner's list 
[38] : (iii) double-confluent Heun equations for the inverse fourth-power potential 
V{r) = V\r^ A + V^r -3 + V%r~ 2 + V^r^ 1 [371 EH], and for the asymmetric double- 
Morse potential given by Zaslavskii and Ulyanov [44j : (iv) confluent Heun equations 
for the trigonometric and hyperbolic potentials given by Ushveridze [39]; (v) general 
Heun equations in the Darboux elliptic form [9] for the first and second Ushveridze's 
elliptic potentials |39|. 

Finally we write a potential given by Cho and Ho [5J, namely, 

V(u) = -\ sinh 2 u - ( I 2 - ]) -4l-. u G (-°°> °°)» (^ = 1, 2, 3, • • • ) (A.7) 
4 \ 4/ cosh u 

where b is a positive real constant. This is a bottomless potential in the sense that 
V{u) — ► —00 when u — > ±00. If b < 4£ 2 — 1, it is an inverted double-well potential; 
if b > 4£ 2 — 1, the potential is similar to the one of an inverted oscillator. For this 
potential the Schrodinger reduces to the Whittaker-Ince limit ([6]) of the confluent 
Heun equation. 

Appendix B. Confluent Hypergeometric Functions 

The regular and irregular confluent hypergeometric functions, denoted respec- 
tively by <5(a, c;y) and ^>(a,c;y), satisfy the Kummer transformations 

<f>{&,c;y) = e y $(c-a,c;-y), ^(a,c;y) = y 1-c *(l + 0- c,2 - c;y). (B.l) 

The behaviour of of ^(a, c; y) and <3?(a, c; y) when y — ► 00 is given by 

i-lvi ( 3tt ^ 3tt\ 



lim*(a,c;y)~y-°[l + 0(|y|- 1 )] ) -— < argy < — (B.2) 



lim 3>(a, c;y) = < 

y— >oo 



' ^\ e y y ^[i + o(\ y \-% (%>0) 
r(o) 

' !ri (- y )-«[l + (|y|-l)], (%<0). 



(B.3) 



, T(c — a) 
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These can be used to get the limits of the series expansions as y —* oo. For $ty = 0, 
the limit of <3?(a, c;y) is a combination of the limits on the right-hand side of the 
previous expression [12J. 

To find the expansions in series of irregular confluent hypergeometric functions 
for the CHE, in addition to Eq. (123D . we need the relations 



yd^pvl = (i _ c)9(flf c; y) + (c _ a _ i)$( 0j c _ i ; y)j 

^yl = ^(a,c;y)-^(a,c+l;y), 

(c-a- l)tf (a, c - 1; y) + (1 - c - y)^(o, c; y) + y^(a, c + 1; y) = 0, 



(B.4) 



where the last relation results from the first and second ones. These relations also 
hold for the functions §(a,b;y) defined in Eq. (|32p . On the other hand, if c — ► oo, 
while a and y remain bounded, we have [12] 



y(a,c;y) = (c)- a (-l)" a + (V^F/r(o)) c ^-(3/2) y i-c e y-c 



(c)" a 



'9tt / r \ c + a -( 3 / 2 ) 
(_!)-« + ^ (_L) a-(l/2) e y+a-(3/2) 

r(o) Vey/ 



[i + o( 



l + O 



which, in conjunction with the last realion of Eq. (1B.4J) . is useful in the study of the 
convergence of the series solutions for the CHE. Notice that this equation implies 



lim ^(a,c;y)~(-c)- a [l + 0(| C |- 1 )], 



lim vf^qy)-^ (±\*" "' ,-"„'<-4,.</+»- "^ 



T(a) \ey 



C -"y""2 e i 



l + O 



(B.5) 
(B.6) 



The Whittaker-Ince limit of the expansions in series of confluent hypergeometric 
functions for the CHE and DCHE yields expansions in series of Bessel functions for 
Eqs. ([6]) and (|7j) by means of the limits [12] 



lim <!> (a,c; — 



r(c)y( 1 - c )/ 2 J c _ 1 (2Vy) ; 



(B.7) 



lim 

a— >oo 



r(a + l -c) * (a,c; 



2y( 1 - c )/ 2 ET c _ 1 (2Vy), 



(Bi 



lim 

a— >oo 



r(a + 1 -c) * (a,c; 



'-ivre i7rc y( 1 - c )/ 2 ^i ) 1 (2Vy),Im y > 0, 



^-^(1-^/2^) 



(B.9) 



^ 1 (2Vy),Imy<0, 



where Ja denotes the Bessel functions of the first kind of order A, if \ are modified 



A 



rf and flf 
Hankel functions [13], respectively. Some relations useful for Sec. 3.2 are 



Bessel functions of the third kind, whereas H[ ' and H\ ' denote the first and second 



H ( y x {x) = e™ x H ( V(x), H (2) x (x) = e~ OTA #f } (x), 

K x (-ix) = \me^ iX H^\x), K x {ix) = -\me~^ iX Hf\ 



(B.10) 
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